Quasi-Monte Carlo and Multilevel Monte Carlo Methods for Computing Posterior Expectations in Elliptic Inverse Problems∗
نویسندگان
چکیده
We are interested in computing the expectation of a functional of a PDE solution under a Bayesian posterior distribution. Using Bayes’s rule, we reduce the problem to estimating the ratio of two related prior expectations. For a model elliptic problem, we provide a full convergence and complexity analysis of the ratio estimator in the case where Monte Carlo, quasi-Monte Carlo, or multilevel Monte Carlo methods are used as estimators for the two prior expectations. We show that the computational complexity of the ratio estimator to achieve a given accuracy is the same as the corresponding complexity of the individual estimators for the numerator and the denominator. We also include numerical simulations, in the context of the model elliptic problem, which demonstrate the effectiveness of the approach.
منابع مشابه
Sequential Monte Carlo methods for Bayesian elliptic inverse problems
In this article we consider a Bayesian inverse problem associated to elliptic partial differential equations (PDEs) in two and three dimensions. This class of inverse problems is important in applications such as hydrology, but the complexity of the link function between unknown field and measurements can make it difficult to draw inference from the associated posterior. We prove that for this ...
متن کاملDomain Decomposition Solution of Elliptic Boundary-Value Problems via Monte Carlo and Quasi-Monte Carlo Methods
Domain decomposition of two-dimensional domains on which boundary-value elliptic problems are formulated, is accomplished by probabilistic (Monte Carlo) as well as by quasi-Monte Carlo methods, generating only few interfacial values and interpolating on them. Continuous approximations for the trace of solution are thus obtained, to be used as boundary data for the sub-problems. The numerical tr...
متن کاملMultilevel higher order Quasi-Monte Carlo Bayesian Estimation
We propose and analyze deterministic multilevel approximations for Bayesian inversion of operator equations with uncertain distributed parameters, subject to additive gaussian measurement data. The algorithms use a multilevel (ML) approach based on deterministic, higher order quasi-Monte Carlo (HoQMC) quadrature for approximating the high-dimensional expectations, which arise in the Bayesian es...
متن کاملLikelihood-informed dimension reduction for nonlinear inverse problems
The intrinsic dimensionality of an inverse problem is affected by prior information, the accuracy and number of observations, and the smoothing properties of the forward operator. From a Bayesian perspective, changes from the prior to the posterior may, in many problems, be confined to a relatively lowdimensional subspace of the parameter space. We present a dimension reduction approach that de...
متن کاملMultilevel Monte Carlo methods and applications to elliptic PDEs with random coefficients
We consider the numerical solution of elliptic partial differential equations with random coefficients. Such problems arise, for example, in uncertainty quantification for groundwater flow. We describe a novel variance reduction technique for the standard Monte Carlo method, called the multilevel Monte Carlo method. The main result is that in certain circumstances the asymptotic cost of solving...
متن کامل